Hindrance of heavy-ion fusion due to nuclear incompressibility 
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We propose a new mechanism to explain the unexpected steep falloff of fusion cross sections 
at energies far below the Coulomb barrier. The saturation properties of nuclear matter are caus- 
ing a hindrance to large overlap of the reacting nuclei and consequently a sensitive change of the 
nuclear potential inside the barrier. We report in this letter a good agreement with the data of 
coupled-channels calculation for the ^*Ni-|-^*Ni combination using the double-folding potential with 
M3Y-Reid effective N — N forces supplemented with a repulsive core that reproduces the nuclear 
incompressibility for total overlap. 
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In recent years a new phenomenon observed in the fu- 
sion of several medium-heavy dinuclear systems Q, 0, 
1^1 ^ 1^ has presented a challenge to the theoretical 
understanding of the reaction mechanism at deep sub- 
barrier energies. From the inspection of capture cross 
sections for various combinations of colliding systems 
such as ssNi-h^SNi, <^^m+ 64Ni, ^ONi-^ ^^Y, '^'^m+^°°Mo, 
902r+90Zr, ^°Zt+^^Y and ^°Zt+^^Zt, an unexpected hin- 
drance was observed at sub-barrier bombarding energies 
below a certain threshold Es that varies from system to 
system. The authors of Ref. remarked that the ex- 
ploration of the hindrance phenomenon is only in its ini- 
tial stage and an identification of the underlying physical 
cause is still missing. The standard theoretical approach 
to treat capture reactions, i.e., the quantum tunneling 
through the relative barrier of the dinuclear system cou- 
pled to different low-lying collective degrees of freedom 
(vibrations and rotations) was unable to explain the steep 
falloff in the cross sections. The coupled-channels (c.c.) 
approach failed to describe the dynamics of the fusion 
process below a certain value of the bombarding energy. 
The source of this phenomenon is not only interesting 
for the understanding of the reaction mechanism but may 
have essential consequences for nuclear reactions that oc- 
cur in stars . It could imply that the synthesis of heavy 
elements is hindered below a certain energy threshold. 

To resolve the enigma several hypotheses have been 
proposed. An issue debated by several authors was the 
large diffuseness needed to fit high-precision fusion data 
[^□l. A suggestion that is very close in spirit to the 
present work was put forward in Q. The authors made 
the observation that potentials such as the Winther- 
Akyiiz provide reliable barriers but the fact that they 
cannot reproduce the data far below the barrier is a sig- 
nature that the ion-ion potential has another form in the 
inner part of the barrier. Using simple arguments (no 



channel coupling, WKB calculation of the transmission 
coefficient) they pointed out that the steep falloff in the 
tunneling probability is related to the disappearance of 
the classically allowed region below a certain energy. If 
this is true, it would mean that we are confronted with 
the existence of a shallow pocket of the potential inside 
the barrier. 

Another interesting observation Q made in relation 
to the data was that the astrophysical 5-factor, as de- 
fined by Burbidgc et al. [loj . develops a maximum for 
the systems enumerated above, which, again, cannot be 
reproduced by c.c. calculations using the deep Winther- 
Akyiiz potential. For the systems mentioned earlier the 
maximum occurs at an energy 20-30 MeV above the com- 
pound nucleus ground state. 

We advocate in this work the idea that in order to an- 
alyze the fusion data the double-folding potentials may 
avoid the systematic failure of other potentials. To reach 
this goal we have to make some amendments to the 
schemes that are usually employed; see for example Ref. 
flll | . First of all we need to take into account the satura- 
tion of the nuclear matter. Second, the neutron and pro- 
ton content of the various dinuclear combinations have to 
be included in the potential, a fact which is often over- 
looked or only indirectly accounted for in the Woods- 
Saxon parametrization. We show that by properly ad- 
dressing these issues light can be shed on the hindrance 
in sub-barrier fusion. 

There are several facts pointing toward the existence of 
shallow pockets in the relative heavy-ion potential used 
to describe various nuclear reactions. The first evidence 
came from the discovery of resonant structures in colli- 
sions of light nuclei, the best known example being the 
sharp peaks in the bombarding energy dependence of the 
excitation curves measured in the ^^C+^^C scattering by 
Bromley et al. |l2l | . These resonances would then resem- 
ble states in a quasi-molecular potential well. 

It became clear in the last two decades that the res- 
onant behavior observed not only in ^^C+^^C but also 
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O, i^C+ "C, 1^0, ^-^Mg, is not 



an isolated phenomenon occurring only in a few a lighter 
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systems. It persists even in heavier colliding; systems, 
such as the '^^Mg+ ^^Mg [l^ and "^^81+ ^^Si Man- 
ifestation of clusterization in relation to quasi-molecular 
pockets is also known for various heavy nuclear systems 
such as cluster radioactivity or cold fission of ac- 
tinides [lij . 

To simulate the appearance of shallow pockets several 
recipes have been proposed in the past : a) A Gaus- 
sian added to the conventional Woods-Saxon potential 
was used in Refs. and [13 in order to fit the re- 

action cross sections observed in ^'^C+ ^^C, ^'^C+ ^^O 
and ^^O; b) potentials computed within the den- 

sity functional method 0, E(l| and c) the well-known 
proximity potential (sec Ref. for the 1977 and jl^l 
for the 2000 versions). 

A repulsive core can be also simulated by folding two 
nuclear density distributions with an effective nucleon- 
nucleon {N — N) interaction. Double folding potentials, 
without repulsive cores, introduced by Satchler and Love 
[i^l, are accurate only in the tail region of the nucleus- 
nucleus potential where the density distributions are only 
gently overlapping and thus the assumption of "frozen 
density" is less questionable. However this assumption 
ignores any readjustment due to the mutual excitations 
of the nuclei or the Pauli exclusion principle for strong 
overlap. Density dependent interactions have been used 
in recent times to simulate the saturation of nuclear mat- 
ter for a+nucleus scattering (see Ref. and references 
therein). For the target-projectile combinations used in 
the experiments with the ATLAS facility at Argonne 
Q, 1^, 1^ Q these double-folding potentials do not nec- 
essarily apply because the resulting pockets are still too 
deep and the barriers too thin below a certain energy. 
For these reasons another approach is considered to in- 
corporate the effect of the nuclear incompressibility. 

To compute the double-folding potential we use the 
Fourier technique as expounded in a previous publication 
|25l |. Since we are interested in performing dynamical 
calculations we consider two heavy ions with one-body 
deformed densities pi and p2, subjected to vibrational 
fluctuations, and center of masses separated by the dis- 
tance R. Then the interaction between these two ions 
can be evaluated as the double-folding integral of these 
densities 

V{R) ^ j dri J dr 2 pi{ri)p2ir2)v{r 12, pi,p2), (1) 

where ri2 = i? + r2 — Ti. In the above formula there 
is allowance for a density dependence of v as discussed 
previously. The central part of the effective N — N poten- 
tial V contains a direct part that also depends on isospin 
since in all the cases of interest N ^ Z, and an exchange 
part that takes into account the effect of antisymmetriza- 
tion under exchange of nucleons between the two nuclei. 
For the density independent part of the effective nucleon- 
nucleon force we use the Michigan-3-Yukawa(M3Y)-Rcid 



parametrization [i l lis l li^. 

The density independent M3Y interactions are cor- 
rectly predicting the potential for peripheral collisions. 
However reactions sensitive to the potential at smaller 
distances are not reproduced j^^l . To cure this deficiency 
the ion-ion potential should also contain a repulsive core. 
For the repulsive part we take a double-folding integral as 
in Eq. but use densities with smaller diffusivity (a < 
0.4) and the N — N interaction is modeled by a zero-range 
form with strength T4orc, following the suggestion of Ref. 
[27j . In what follows we refer to the used potential (in- 
cluding the Coulomb part that is also calculated via the 
double-folding procedure) as the M3Y-|-repulsion. 

To derive the properties of the short-range repulsive 
core, we note that an overlapping region with doubled 
nucleon density is formed once the distance R between 
the nuclei becomes less than Rp+Rt, where Rp and Rt are 
the nuclear radii along the collision axis. The doubling 
of the density increases the energy of the nucleons in the 
overlapping region. In the case of complete overlap (for 
i? = 0) the increase of the interaction energy per nucleon 
is, up to quadratic terms in the normal density, 

where in the last equality we use the proportionality of 
the incompressibility K of normal nuclear matter to the 
curvature of the energy per nucleon |2^ . Since K is usu- 
ally not measured directly but deduced from isoscalar 
giant monopole or dipole resonances, and since there are 
conflicting results coming from the cross sections of the 
corresponding experimental data (see and references 
therein) we use instead of a universal value the predic- 
tions of the Thomas-Fermi model |33 as a function of the 
relative neutron excess S = (/?„ — Pp)/p of the compound 
(fused) nucleus. Eventually the strength of the repulsive 
core, V^orc, is determined by assuming that AV must 
be identified with the value of the heavy-ion interaction 
potential at the coordinate origin R = 0. 

The parameters of the Fermi-Dirac density distribution 
entering in the double-folding potential ^ are calculated 
as follows: For the proton density distribution we use the 
parameters resultin g fro m a compilation of elastic elec- 
tron scattering data [31|, i.e. rop=1.065 fm and ap=0.575 
fm. For the neutron distribution we took the liberty 
to vary the parameters under reasonable limits, taking 
into account the following constraints: a) we deal with 
a moderately neutron rich nucleus and b) the barrier of 
the M3Y-|-repulsion potential should be close to the one 
predicted by the Winther-Akyiiz potential, since this po- 
tential gives a good description of the data in the barrier 
and high-energy region Thence we take ro„=1.085 
fm and a„=0.534 fm. To obtain the incompressibility of 
A'=228 MeV, predicted by the Thomas-Fermi model for 
the compound nucleus ^^^Ba [s^l, we use a strength of 
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FIG. 1: Various spherical ion-ion potentials for '"'Ni+^Ni. 
The solid curve is the potential employed in the present work. 
The curve with small dashes is the Winther-Akyiiz potential 
used in Q, |^. The dashed strip corresponds to exper- 

imental boundaries of the threshold energy Eg. 

14oro=496 MeV, and a difFuseness of the repulsive part 
of the density distribution of a wO.4 fm for both protons 
and neutrons. 

To illustrate the shape of this potential for a case of 
interest we compare in Fig. ^the spherical heavy-ion po- 
tential for the symmetric dinuclear system ^"'Ni+^^Ni to 
different potentials that have been used in the past to ex- 
plain fusion da ta ( Winther-Akyiiz 0]) or fusion barriers 
(Proximity '77 |2l|). The thickness of the dashed region 
in Fig. ^reflects the uncertainty in the energy Eg, where 
the experimental iS-factor has a maximum. 

As a study case we consider the system ^^Ni-l-^^'Ni 
and treat the fusion dynamics within the c.c. method 
with linear and quadratic couplings in the quadrupole 
and octupole vibrational amplitudes |32l |. The c.c. equa- 
tions, which include couplings to the low-lying 2+ and 3~ 
states, mutual excitations, and two-phonon quadrupole 
excitations, are solved in the rotating frame approxima- 
tion with the usual outgoing-wave boundary conditions 
at large distances. Ingoing-wave boundary conditions are 
imposed at the radial separation where the heavy-ion po- 
tential develops a local minimum, Vi„in ~ 85 MeV. As a 
result of this scheme, the calculated cross section will van- 
ish for center-of-mass energies E < Vmin- Naturally one 
would expect a non-zero cross section even below Vmin, 
since the compound nucleus can in principle be formed 
down to a threshold energy of E fa 49 MeV. 

In Fig. 121 the experimental excitation function of the 
fusion reaction ^^Ni -I- ^''Ni ^^^^Ba is compared with 
the results using the Winther-Akyiiz potential as in Ref. 
(dotted line) and with the M3Y-|-repulsion potential 
(solid line), in both cases using the same recipe for the 
c.c. calculations. The cross section for no couplings is 



FIG. 2: Experimental fusion excitation function for the sys- 
tem ^''Ni-|-^'*Ni 1^ is compared to various calculations de- 
scribed in the text. 



also shown (dashed curve) when the Winther-Akyiiz po- 
tential is used. It is seen in Fig. |21 that the agreement 
with the data, when using the M3Y+repulsion potential, 
is much better than the one provided by the Winther- 
Akyiiz potential. The " M3Y+ repulsion" excitation func- 
tion has the right shape, not only because the potential 
attains a higher-lying minimum but also because the cur- 
vature of the barrier is different, producing a thicker clas- 
sically forbidden region. Taking all the experimental data 
points into account we obtain a minimum P^r point 
of x^/N = 10.1 for an overall energy shift A£;=0.9 MeV 
of the calculated cross sections when the Winther-Akyiiz 
potential is used. For the M3Y-|-repulsion we get a much 
better result with a minimum /N = 0.86 for an energy 
shift of A£;=0.16 MeV. 

The S'-factor representation of the ^"'Ni-l-^^Ni data 
(solid points) is compared in Fig. |31with the two c.c. cal- 
culations, based on the M3Y+repulsion potential (solid 
curve) and the Winther-Akyiiz potential (dashed curve) . 
The clear maximum in the S'-factor is reproduced only by 
the M3Y+repulsion potential. At this point one should 
recall the experience gained in the sixties on molecular 
resonances. As shown in [s^l the S-f actor exhibits a se- 
quence of quasi-molecular resonances for lighter systems. 
In the present case we get a maximum that is too broad to 
be assigned to a resonance, the curvature in the S'-factor 
being explained by the shallow pocket in the potential. 
The maximum of the theoretical curve corresponds ap- 
proximately to the maximum of the data points. 

A similar study was performed on another case, 
^*Ni-|-^*Ni, for which data are available from an older 
experiment with the smallest cross sections in the range 
of mb [31 ■ Even in this case we confirm the trend of 
the excitation functions, namely a faster decrease than 
expected due to the existence of a repulsive core. 
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FIG. 3: The experimental S-factor of ^''Ni+'^'»Ni 0] (solid 
points) is compared to the coupled-channels calculations 
performed with the M3Y+repulsion (solid curve) and the 
Winther-Akyiiz (dashed curve) potential. 

Thus, the understanding of the experimental data, at 
least for the case ^^Ni+^^Ni, requires a modified shape 
of the potential inside the barrier. Both a thicker bar- 
rier and a shallower pocket are needed to explain the 
measured cross sections and these features are naturally 
explained by the M3Y+repulsion heavy-ion potential. If 
the actual potential pocket would be deeper than the po- 
tential used by us, e.g. the double-folding with density 
dependent M3Y interactions , then one would expect 
within the energy range between the barrier top and the 
threshold energy Es the appearance of a resonant struc- 
ture of the form encountered for lighter quasi-molecular 
systems 113 . Apparently there is no indication of reso- 
nant structures in the excitation functions of ^^Ni-|-^'*Ni 
or in the other systems studied in Refs. flHSSIS- 
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